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Abstract 

The state-dependent diffusion, which concerns the Brownian motion of a particle in inhomoge- 
neous media has been described phenomenologically in a number of ways. Based on a system- 
reservoir nonlinear coupling model we present a microscopic approach to quantum state-dependent 
diffusion and multiplicative noise in terms of a quantum Markovian Langevin description and an 
associated Fokker-Planck equation in position space in the overdamped limit. We examine the 
thermodynamic consistency and explore the possibility of observing a quantum current, a generic 
quantum effect, as a consequence of this state-dependent diffusion similar to one proposed by 
Biittiker [Z. Phys. B 68, 161 (1987)] in a classical context several years ago. 
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I. INTRODUCTION 



Almost three decades ago Landauer^ explored the problem of characterizing nonequilib- 
rium steady states in the transition kinetics between the two locally stable states in bistable 
systems. His main idea was that the relative stability of a particle diffusing in a bistable 
potential can be altered by an intervening hot layer which has the effect of pumping par- 
ticles from a globally stable region to a metastable region. No detailed consideration of 
immediate neighbourhood of the two states is important. In formulating the problem in 
terms of diffusion equation it was realized that one needs state dependence of diffusion for 
a correct description of the effect and more generally a careful analysis of the problem of 
diffusion in inhomogeneous media in a wider context was necessary. This was carried out 
by van Kampen 3 and others^ in eighties. An important consequence of state-dependent 
diffusion or noise as suggested by Biittiker— is the generation of current, in absence of any 
externally applied fields, which occurs in presence of periodic diffusion of a particle in a 
spatially periodic potential with same periodicity but differing in phase. This rectification 
of state dependent noise resulting in a directed transport and state-dependent diffusion play 
important role in several areas of condensed matter physics on the mesoscopic scal o 5 i 6 i 7 i 8 
and furthermore in ratchet problem o 9 ^ ! 11 ^ 2 in a wider perspective. 

The physics of state-dependent diffusion can be described phenomenologically in a num- 
ber of ways. As reported in the literature^ 3 - the diffusion term for Brownian particle may 
assume several forms, notably -^D(q)-^P(q,t) or -^D(q)P(q 7 t) or the other like -§^D 9P ^ 
supplemented by 'thermal potential' or state-dependent drift term. Here D(q) refers to dif- 
fusion coefficient and P(q,t) is the probability distribution function for the particle. Two 
important points are now noteworthy. First, the phenomeno logical forms of diffusion co- 
efficient being different, it is easy to realize that they do not have a common microscopic 
Hamiltonian origin. Thus the physics of diffusion in inhomogeneous media is somewhat 
model-dependent^ 3 " and the search for a "correct" form of diffusion remains a debatable 
issue. Second, the diverse forms notwithstanding, the generalization of Boltzmann factor, 
exp(— V(q)/kBT) (which governs the system at thermal equilibrium) for state-dependent 
diffusion in the steady state assumes a common structure, 



P st (q) ~ exp[-0(g)] 
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with 0(g) = f{f V d \^ dq\ V(q) being the potential field. The steady state distribution 
implies that the effective potential 0(g) is nonlocal in space. The generality in the 
structure of 0(g) is such that it may include the spatial variation of temperature, diffusion 
or drift coefficient as specific CclSCS clS considered separately by several autho ra 3 ' 4 ' 13 . In the 
Langevin scheme of description, on the other hand state-dependent diffusion has received 
attention under "multiplicative noise" 1 ^i 1 7i 1 ?i 1 ?i 2 9i 21 i 22 i 2 ?i 24 . The microscopic origin of multi- 
plicative noise within the framework of standard paradigm of system-reservoir Hamiltonian 
that includes a variety of model calculations is the nonlinear coupling between the system 
and bath coordinates which leads to nonlinear dissipation. A thermodynamically consistent 
approach in this context was put forward in early eighties by Lindenberg and coworkers^. An 
exact Fokker-Planck equation for time and space dependent friction was derived by Pollak 
et aim several years ago. Along with these formal development a 13 ! 14 ! 1 ^ 1 ^ 1 ^ 1 ^ 1 ^ 2 ^ 21 ! 22 ! 23 ' 24 , 
the theories of multiplicative noise have found wide applications in several areas, e.g., acti- 
vated processes* 2 ^ stochastic resonance^, laser and optics^, signal processing 28 , fluctuation- 
induced transpor t 29 ' 30 , noise-induced transitions* 3 ^ etc. 

In this paper we address the problem of Langevin equation with multiplicative noise 
and state-dependent diffusion for a thermodynamically closed system in a quantum me- 
chanical context. Although the quantum-mechanical system- reservoir linear coupling model 
for microscopic description of additive noise and linear dissipation which are related by 
fluctuation-dissipation relation is well-known over many decades in several field o 32 ' 33 , the 
nature of nonlinear coupling and its consequences have been explored with renewed interest 
only recently. For example, it has been observed that the quantum dissipation can reduce 
the appearance of metastable state and barrier drift in a double well potential 34 . Tanimura 
and co-workers*^ have extensively used nonlinear coupling in modeling elastic and inelastic 
relaxation mechanisms and their interplay in vibrational and Raman spectroscopy. The role 
of inhomogeneous dissipation in reducing quantum decay rate has also been explored very 
recently 3 ^. Based on a coherent state representation of quantum noise operator and Wigner 
canonical thermal distribution for harmonic bath oscillators 37 nonlinearly coupled to a sys- 
tem we would like to develop a microscopic approach to quantum state-dependent diffusion 
and quantum multiplicative noise. Specifically, our aim is: 

(i) to derive a quantum Smoluchowski equation for state-dependent diffusion on the basis 
of a system-reservoir model and analyze the role of nonlinear coupling in state-dependence 
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of quantum diffusion and dissipation. 

(ii) to seek for a correspondence between this quantum equation and its phenomenological 
counterpart. 

(iii) to derive a quantum generalization of Boltzmann factor for state-dependent diffusion 
in the steady state and check its thermodynamic consistency. 

(iv) to explore the possibility of observing a directed transport as a consequence of state 
dependent quantum diffusion in the spirit of Biittiker as an immediate application. 

The outlay of the paper is as follows: In Sec. II we develop the scheme of quantum 
Brownian motion for multiplicative noise on the basis of a nonlinearly coupled system- 
reservoir model of Zwanzig-type form within a Markovian description. Sec. Ill is devoted 
to quantum Smoluchowski equation for state-dependent diffusion which corresponds to a 
specific form of a classical phenomenological equation. We derive a quantum mechanical 
generalization of Boltzmann factor for state-dependent diffusion. In Sec. IV we carry out 
an application to a system with spatially periodic potential and periodic diffusion function 
with same periodicity but differing in phase to demonstrate a nonzero current. The paper 
is concluded in Sec.V. 

II. QUANTUM MULTIPLICATIVE NOISE 
A. General aspects; Langevin equation 

We consider a particle of unit mass coupled to a medium comprised of a set of har- 
monic oscillators with frequency 0Uj. This is described by the following system-bath 
Hamiltonia n 34 i 35 i 36 . 



Here q and p are the coordinate and momentum operators of the particle and the {xj,pj} 
are the set of coordinate and momentum operators for the bath oscillators with unit mass. 
The system particle is coupled to the bath oscillators nonlinearly through the general cou- 
pling terms ^-f{q)- Cj is the coupling strength. The Hamiltonian Eq. (|2.1|) is different from 
Zwanzig 38 form of system-bath Hamiltonian where the coupling is linear with respect to 
system coordinate. The classical counterpart 14 of the form (|2.1|) is known for more than 
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two decades and also the nonlinear coupling in quantum system has been studied in several 
context o 34 ! 36 . The potential V(q) is due to external force field for the system particle. The 
coordinate and momentum operators follow the usual commutation relations [q,p] = ih and 
[%,Pfc] = ihdjk- The presence of counter term in the Hamiltonian ensures that the potential 
V(q) felt by the particle does not get modified due to heat bath. 

We now use Eq. (|2.1|) to obtain the following dynamical equations for the position and 
momentum operators: 

• 9H . . 

q=^r=P (2.2) 



where the dot(.) indicates derivative with respect to time and the prime (/) refers to 
derivative with respect to q. 

Similarly we have the dynamical equations of motion for the bath oscillators (j = 1, 2, 3...) 

dH 

Pj = ~dx~ = ~ + ° jf ^ ^'^ 
To eliminate the bath degrees of freedom from the equations of motion of the system we 
first obtain a solution for the position operator Xj by formally solving the above equations 
and then make use of the solution in Eq. (j2.3|) followed by some rearrangement. This yields 
the generalized operator Langevin equation for the system particle. 



i(t) = p{t) (2.6) 

m = -vim) - fw)) f fm')Ht - t'W)dt' + rmmt) (2.7) 



where the noise operator fj(t) and the memory kernel j(t) are given by 



" f 2 ""x 2 

A L \. J } J J 



(2-8) 
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and 



7(t) = ^2^2 cosojjt (2.9) 

It is clear from the operator Langevin equation Eq. (j2.7j) for the system that the noise 
operator is multiplicative and the dissipative term is nonlinear with respect to system co- 
ordinate due to the nonlinear coupling term in the system-bath Hamiltonian. In the case 
of linear coupling, i.e., f(q) — q the Eq. (j2.7j) reduces to a quantum generalized Langevin 
equation 32 in which the noise term is additive and the dissipative term is linear. 

Since the system is thermo dynamically closed, i.e., the fluctuation and the dissipation 
originates from the same origin, the detailed balance condition must be satisfied. The noise 
properties of fj(t) can be derived by using suitable canonical thermal distribution of bath 
coordinates and momenta operators at t — to obtain; 

m))QS = (2.10) 



+ W)v(t))Qs = \H% h ( coth 2Ft) cosu ^ - £ ') ( 2 - n ) 
j j \ b 

Here (....) qs implies quantum statistical average on bath degrees of freedom and is defined 

as 

{ Q )QS = Trd ^(-H bath /k B T) (212) 
Trexp(-H bath /k B T) 



bath 



for any bath operator O {\^-Xj — f(q)}, {pj} j , where H t 

\f j2 c 1 

J2j ~2 + ~ ^-f{q)) 2 a ^ t = 0. By trace we mean the usual quantum statis- 

tical average. Eq. (j2.11j) is the fluctuation-dissipation relation expressed in terms of noise 
operators appropriately ordered in the quantum mechanical sense. 

In the Markovian limit the generalized quantum Langevin equation Eq. (j2.7|) reduces to 
the form 



q(t) = p(t) (2.13a) 

m = -v(q(t)) - r irm))? m + nmw) (2.13b) 
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where T is dissipation constant in the Markovian limit. 

To construct a c-number quantum Langevin equation we procee d 39 ' 40 ! 41 ! 42 ' 43 as follows. 
We carry out a quantum mechanical average of Eq. (j2.13aj) and Eq. (j2.13b|) to get 

q = p (2.14) 
V = -{V'(q)) - T([f'(q)} 2 p) + (f'(q)f)(t)) (2-15) 

where q = (q) and p = (p). The quantum mechanical average (...) is taken over the 
initial product separable quantum states of the particle and the bath oscillators at t = 
0, |0){|q!i)|q!2)....|q;jv)}. Here \<p) denotes any arbitrary initial state of the system and 
\aj) corresponds to the initial coherent state of the j-th bath oscillator. \aj) is given by 
lay) = exp(— |a||/2) Yl^=o{ a ^ 3 1 \fn?)\ n j) i a j being expressed in terms of the mean values 

to! 2 

of the shifted coordinate and momentum of the j-th oscillator, {-^-(^(O)) — (/(<?(0)))} = 
J -^-.{oij + a*) and (pj(0)) = \J ^-(a* — aj), respectively. 

Since the quantum mechanical average is taken over the initial product separable quantum 
states of the particle and the bath oscillators the Eq. ()2. 14112. 15J) can be written as: 

q = p (2.16a) 
P = -WW - r([/'(g)] 2 p) + (f'(q)Ht) (2.16b) 

where r](t) = (fj(t)), rj(t) is now a classical-like noise term, which, in general, is a non-zero 
number because of the quantum mechanical averaging and is given by 

r r to 2 ^ c 2 

V(t) = E - @i(°)> - (/W(°))> 4 C0S ^+ T (Pi(0)>sin^ (2.17) 
j L I j J u j u j J 

To realize rj(t) as an effective c-number noise we now introduce the ansat a 37 i 39 i 40 i 41 i 42 i 43 

U) 2 

that the momentum (ftj(0)) and the shifted coordinates {-^- (xj(0)) — (/(<?(0)))} of the bath 
oscillators are distributed according to a canonical distribution of Gaussian form as: 



^•({f (%(o)) - </(g(o))», fo(o)) 



2 2 

[fe(o)> 2 +^{^ fe(0)M/(g(0))}} 2 ] 
A^exp { i—L—— \ (2.18) 

1 > 2/^(n J K)+i) i 
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so that for any quantum mechanical mean value, Oj f{^r (xj(Q)) — (/(<?(0)))}, (Pj(0))j 
of the bath operator, its statistical average (...)s is 



{0 3 ) s = / Oft d<Pi(0)) d{-^(0)> - (/(g(0)))} (2.19) 

Here nj{ujj) indicates the average thermal photon number of the j-th oscillator at the 
temperature T and nj{ojj) = [exp(^— |?) — 1] _1 and A/" is the normalization constant. 

The distribution Vj (Eq. (|2.18|l ) and the definition of statistical average Eq. (j2.19j) imply 
that c- number noise r](t) must satisfy 

(r)(t)) s = (2.20) 

(V(t)v(t'))s = \H% (coth ^) coautf - (2.21) 

which are equivalent to (|2.1U|) and ()2.11j) . respectively. 
In the Markovian limit the noise correlation becomes 



(V(t)v(t'))s = 2D 5{t-t') (2.22a) 
D = hhuo (u(lj ) + ij (2.22b) 

where c^o is the average bath frequency and the spectral density function is considered in 
the Ohmic limit. 

The Eqs. fl2.20H2.21l) imply that the c-number noise r](t) is such that it is zero centered and 
satisfies the standard fluctuation-dissipation relation as expressed in Eq. fj2.llj ) . It is easy to 
recognize that the ansatz (|2.18j) is a canonical thermal Wigner distribution function for a 
shifted harmonic oscillator (obtained as an exact solution of Wigner equation^ for harmonic 
oscillator) which always remains a positive definite function. A special advantage of using 
this distribution function is that it remains valid as a pure state nonsingular distribution 
function even at T = 0. At the same time the distribution of quantum mechanical mean 
values of the bath oscillators reduces to classical Maxwell-Boltzmann distribution in the 
thermal limit fuv <C kgT. Furthermore this procedure allows us to bypass operator ordering 
prescription of Eq. (|2.11j) for deriving noise properties of the bath oscillators and to identify 
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r)(t) as a classical looking noise with quantum mechanical content. We also mention that 
instead of Wigner function it is also possible to employ Glauber-Sudarshan distribution 
function to derive quantum fluctuation-dissipation relatio n 39 i 40 i 41 i 42 i 43 . 

We now return to Ea. (l2.1baH2T6bj) to add V'(q), T[f'(q)fp and f'(q)i](t) on the both 
sides of Eq. (|2.16b|) and rearrange it to obtain 

q = p (2.23a) 
V = -V'(q) + Qv~ T[f'(q)] 2 p + Q 1 + f'(q) V (t) + Q 2 (2.23b) 

Where, 



Qv = V\q) - (V'(q)) (2.24a) 
Qi = r[{f'(q)] 2 p-(inq)] 2 p)} (2.24b) 
Q2 = V(t) [</'«)> - /'(?)] (2.24c) 

Here Qv represents quantum correction due to nonlinearity of the system potential. 
Qi and Q2 represent quantum corrections due to nonlinearity of the system-bath coupling 
function. This implies that the quantum Langevin equation is governed by a c-number 
noise rj(t) originating from the heat bath characterized by the properties (j2.2()H2.2Tj) and 
the quantum correction terms Qv, Qi and Q 2 are characteristic of the nonlinearity of the 
potential and the coupling function. 

Referring to the quantum nature of the system in the Heisenberg picture we now write 
the system operators q and p as 



q = q + 5q (2.25a) 
p = p + 5p (2.25b) 

where q(= (q)) and p(= (p)) are the quantum mechanical mean values and Sq and Sp 
are the operators and they are quantum fluctuations around their respective mean values. 
By construction (Sq) = (Sp) = and they also follow the usual commutation relation 
[Sq, Sp] = ih. Using ()2.25a|) and (|2.25b|) in V'(q), [f\q)] 2 p and f'(q) and a Taylor series 
expansion in Sq around q it is possible to express Qv, Qi and Q 2 respectively a a 44 ' 45 
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where, 



Qv = -£V n+1 (<zX<^> (2.26) 



n>2 



Q 1 = -T[2pf'(q)Q f +pQ 3 + 2f'(q)Q 4 + Q 5 ] (2.27) 
Q 2 = v{t) Qf (2.28) 



n>2 

m>l n>l 

Q 4 = V^r +1 (g)(5ffi) (2.31) 



n>l 



1 1 



= EEn^^^^x™) ( 2 - 32 ) 



m>l n>l 



Using ()2.27|) and ()2.28|) in Eq. (|2.23b|) we have the c- number quantum Langevin equation 
in the Markovian limit 

q = p (2.33a) 



V = -V'(q) + Q v - r[f(q)} 2 p -2Tp f'(q)Q f -TpQ 3 - 2Tf(q)Q 4 -TQ 5 

+ f'(q)v(t) + QfV(t) (2.33b) 

The quantum Langevin equation is characterized by a classical force term, V, as well as its 
correction Qv- The terms containing T are nonlinear dissipative terms where Qf, Q 3 , Q4 and 
Q5 are due to associated quantum contribution in addition to classical nonlinear dissipative 
term T[f'(q)] 2 p. The last term in the above equation refers to a quantum multiplicative 
noise term in addition to the usual classical contribution f'{q)i](t). It is therefore easy 
to recognize the classical limit of the above equation derived earlier by Lindenberg and 
Seshadri^A. Furthermore quantum dispersions due to potential and coupling terms in the 
Hamiltonian are entangled with nonlinearity. The quantum noise due heat bath on the other 
hand is expressed in terms of the fluctuation-dissipation relation. 
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B. Quantum correction equations 



The quantum correction terms due to nonlinearity of the potential and the coupling 
function in Eq. (j2.33b|) are taken care of to all orders, in principle. In order to calculate 
them and the associated equations explicitly we write the system operators q, p and the 
noise operator f) in the Heisenberg picture as 



q = q + 5q , 
p = p + 5p , 
fj = 7] + 5rj , 



(2.34a) 
(2.34b) 
(2.34c) 



where 5fj is the fluctuation in noise around the mean value rj and also (5fj) = 0. 
We then use the above set of equations in operator Langevin equation in the Markovian 
limit fEq. (J2.13aj) and Eq. (j2.13b|) ) and subtract the c- number quantum Langevin equation 
fEq. ()2.33a|) and Eq. (j2.33b|) ) from the resultant to obtain 



5q = 5p 

8p = -V"{q)8q - ~y n+ \<l) W - (<*<f >] 



^2.35) 



n>2 



m>l n>l 



P 



- r 



2f(q)f(q)Sq + 2f'(q) V -J n+ \q) [5q n - (5q n )} 

n>2 

+ E A V m+1 (^ n+1 ^) [^ m ^ n - w>] 

z — ' ^— ' mini 

n>l 

[f(q)] 2 5p + 2f'(q) ^J n+1 ^) ^ " 

n>l 

+ V V; -lI/-+l( g )/«+l(g) [5<T5q n 5p - (5q m 5q n 5p)\ 
z — ' z — ' m! n! 

m>l n>l 

n>2 

s, i /'(9)+EV +1 (^ n 



n>l 



(2.36) 



We now carry out quantum mechanical average of Eq. ()2.36|) over initial product separable 
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bath states n°^ 1 {|a :) (0))} to get rid of 8fj term. With this the correction equations result 



m. 



5q = 5p 

S'P = -V"(q)8q - E V" +1 ^) " 



n>2 



- r 



2f'(q)f"(q)5q + 2f(q) V -J n+ \q) [5q n - (5q n )] 

^— ' n! 



n>2 



+ E E A V m+1 ^ n+1 (<?) t^ m ^ n - (^ m ^ n )i 

z — ' m! n! 

m>l n>l 



n>l 



1 1 



+ EE ^ v m+1 (^)r +1 (<?) - 

L — ' z — 4 ml n\ 



m>l n>l 



+ r](t) 



n>2 



(2.37a) 



(2.37b) 



The operator equations ()2 . 3 7al and 12 . 3 7b|) are the basis of calculation of quantum correc- 
tion terms Qv, Qf, Q3, Q± and Q 5 to an arbitrary order. More specifically, we are to set up 
the equations for (Sq 2 ), (5q5p + 5p5q), (5p 2 ) in the second order and similarly for third order 
and so on. They are, in general, coupled and the infinite set of hierarchy of equations must 
be truncated after the desired order for practical purpose. The procedure is exactly similar 
to what had been done earlier in the case of additive nois o 39 i 40 i 41 i 42 i 43 . Thus upto third order 
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we may construct the following set of equations from Eq. (|2.37al and I2.37a|) ; 



£(5q 2 ) = (SqSp + 5p5q) (2.38a) 

(JjL 

±{6q5p + 5p6q) = -2 X (q,p)(5q 2 ) + 2(5q 2 ) -T[f'(q)] 2 (5q5p + Sp5q) 

~ C(q,p)(Sq 3 )-2Tf'(q)f(q)(5q 2 5p + 5p5q 2 ) (2.38b) 

j t (5p 2 ) = -2r[f'(q)] 2 (5p 2 )-x(q,p)mp + 5p5q) 

- lc(q,p)(Sq 2 6p + 5p5q 2 ) - 2Tf(q)f(q)(5q5p 2 + 5p 2 5q) (2.38c) 

£{5q 3 ) = -(5q 2 5p + SpSq 2 ) (2.38d) 

at 2 

^(Sp 3 ) = -3T[f\q)] 2 (5p 3 ) - ^ X (q,p)(Sq5p 2 + Sp 2 5q) (2.38e) 

j t (6q 2 6p + 6p5q 2 ) = -2x(q,p)(Sq 3 ) + 2(5qSp 2 + 6p 2 6q) 



T[f'(q)} 2 (5q 2 5p + 5pSq 2 ) (2.38f) 

4x(q,p)(dq 2 $i 

2T[f'(q)} 2 (5q5p 2 + 5p 2 5q) (2.38g) 



-^.{SqSp 2 + Sp 2 5q) = 2(5p 3 ) - 4:x(q,p)(5q 2 5p + 5p5q 2 ) 



where 

X (q,p) = V"{q) + 2Yp f(q)f"(q) - v(t)f"(q) (2.38h) 
C(q,p) = V"\q) + 2Tp f'{q)f" (q) + 2T p [f"{q)} 2 - V (t)f'"(q) (2.38i) 



C. Calculation of quantum statistical averages 

Summarizing the discussions of the last two sections A and B we now see that the quan- 
tum Langevin dynamics can be calculated for a stochastic process by solving the coupled 
Eqs. (2.33a, 2.33b) and (2.38a-2.38g) for quantum mechanical mean values simultaneously 
with quantum correction equations which describe quantum dispersion around these mean 
values. In principle, the equations for quantum corrections constitute an infinite set of 
hierarchy which must be truncated after a desired order, in practice, to make the system 
of equations closed. Care must be taken to distinguish three averages, the quantum me- 
chanical mean (0)(= O), statistical average over quantum mechanical mean (O)s and the 
usual quantum statistical average (O)qs as discussed in Sec. A. To illustrate the relation 
among them let us calculate, for example, the quantum statistical averages (q)qs, (q 2 )QS 
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and {q 2 p)Q S - By (I2.34a|) and ()2.34bj) we write 



q = q + 5q 



(q)qs = (q + Sq)QS 



(2.39) 



(q)s+((Sq))s = (q)s 



(2.40) 



Again 



(<f>QS = ((q + Sq) 2 ) QS 

= (q 2 )s + ((5q 2 ))s 



(2.41) 



In the case of harmonic potential and linear coupling, (Sq 2 ) is independent of q or p so 
that one may simplify (j2.41|) further as 



In Ref.— the explicit exact expressions for (q 2 )Qs and (p 2 )qs have been derived for har- 
monic oscillator in linear system-bath coupling and they are found to be in exact agreement 
with those of Grabert et a/ 46 . For anharmonic potential, however, one must have to use 
(I2.41j) to carry out further the statistical average over (5q 2 ) i.e ((Sq 2 ))s, since (Sq 2 ) is a 
function of stochastic variables q and p according to quantum correction equations. Fur- 
thermore we consider (q 2 p)Qs 



The essential element of the present approach is thus expressing the quantum statistical 
average as the sum of statistical averages of set of functions of quantum mechanical mean 
values and dispersions. Langevin dynamics being coupled to quantum correction equations, 
the quantum mechanical mean values as well as the dispersions are computed simultaneously 
for each realization of the stochastic "path". A statistical average implies the averaging over 
many such "paths" (typically three to five thousands) similar to what is done to calculate 




(2.42) 



(q 2 P)QS = ((q + Sq) 2 {p + Sp)) QS 

= (q 2 p)s + (P (5q 2 ))s + ((5q 2 Sp))s + 2(q(5q Sp)) s 



(2.43) 
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statistical averaging by solving classical Langevin equation. Before leaving this section we 
mention a few pertinent points. 

First, the distinction between the ensemble averaging by the present procedure and by 
the standard approach using Wigner function is now clear. From Eq. (J2.43j) we note that, 
for example, 



where W(q,p) is the Wigner function for the system. (This is not be confused with the 
Wigner function we introduced in Eq. (J2.18|) for the bath oscillators). 

Second, our formulation of the Langevin equation coupled to quantum correction equa- 
tions belongs to quantum stochastic process driven by c-number noise, which is classical-like 
in form. Its numerical solutions can be obtaine d 41 ^ 42 in the same way as one proceeds in a 
classical theory. 

Third, quantum nature of the dynamics appears in two different ways. The heat bath 
is quantum mechanical in character whose noise properties are expressed through quantum 
fluctuation-dissipation relation. The nonlinearity of the system potential and coupling, on 
the other hand, give rise to quantum correction terms. Thus the classical Langevin equation 
can be easily recovered (i) in the limit fkv <C k^T to be applied in the Eq. (j2.21|) so that 
one obtains the classical fluctuation-dissipation relation and (ii) if the quantum dispersion 
terms vanish. 

III. THE OVERDAMPED LIMIT AND THE STATIONARY DISTRIBUTION 

A. The Langevin equation with multiplicative noise under strong friction 

In the case of large dissipation, one eliminates the fast variables adiabatically to get a 
simpler description of the system which is valid in a much slower time scale. This adiabatic 
elimination of fast variables is basically a zeroth order approximation. In this zeroth order 
approximation the number of system variables get reduced. When the Brownian particles 
move in a bath with constant large dissipation this adiabatic elimination of fast variables 
leads to the correct description of the system. This wellknown approximation, known as 
Smoluchowski approximation, results in correct equilibrium distribution. However in pres- 




(2.44) 
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ence of hydrodynamic interaction, i.e., when the fluctuation is position/state dependent 
or equivalently when the noise is multiplicative with respect to system variables the con- 
ventional adiabatic reduction of fast variables does not give the correct description of the 
system. Several years ago Sancho et aU^ had proposed an alternative approach to get a 
correct Langevin equation in the case multiplicative noise system. Based on the Langevin 
equation they carried out a systematic expansion of the relevant variables in powers of 
neglecting terms smaller than 0(T -1 ). Then by ordinary Stratonovich interpretation it is 
possible to obtain the correct Langevin equation corresponding to a Fokker-Planck equation 
in position space. This description leads to the correct stationary probability distribution 
of the system with position dependent friction. 

In order to get the quantum Langevin equation in the overdamped limit we follow this 
procedure. In what follows we discard the quantum correction terms Q4 and Q5 from 
Eq. (|2.33b|) since they involve quantum dispersions Sp and therefore decays exponentially 
in the large damping limit. In this limit these transient correction terms do not affect the 
dynamics of the position which varies in a much slower time scale. So the quantum Langevin 
equation Eq. ()2.33a|) and Eq. ()2.33bj) can be written as, respectively, 



q = p (3.1) 
p = -V'(q) + Q V -T h(q)p + g(q) V (t) (3.2) 



where 



%) = [f'(q)} 2 + V'(q)Qf + Q3 

g{q) = m + Qf 



(3.3) 
(3.4) 



The method of Sancho et al^ is followed further to obtain the Fokker-Planck equation 
in position space corresponding Langevin equation Eq. ()3.2jl 



dP(q,t) 

dt 



d_ 

dq 



V'(q) - Q\ 



+ D 4 
oq 



Th(q) 
' g(q) d g(q) 



[T (h(q))< 



g'(q)g(q) 



r %) d q r %)_ 



P(q,t) 



(3.5) 



16 



In the ordinary Stratonovich description the Langevin equation corresponding to the 
Fokker-Planck Eq. (j3.5|) is given by 



Q=- 



V'(g) - Qv Du g'(g)g(q) | g(g) ^ 



(3.6) 



T h{q) u T {h{q)Y V h(q) 

Eq. (j3.6j) is c-number quantum Langevin equation for multiplicative noise with position 
dependent friction in the overdamped limit (i.e., corrected upto 0(1 /V)). 

In the classical limit, i.e., hj Q <C k B T, h(q) = \f'(q)} 2 , g(q) = f'(q), Qv — and 
Do = T ksT so the c-number quantum Langevin equation reduces to 

1 



Q 



r [f'(q)f 



-V'(q) ~ r k B T + f(q) V (t) 



(3.7) 



which is exactly the form derived by Sancho et al 17 . 

Our next task is to establish the quantum corrections in the overdamped limit. To this 
end we return to Eq. (J2.37bj) . neglect the Sp term and keep the leading order correction 
equation (since in the overdamped limit higher order quantum contributions are small) to 
obtain 



J t 6q= r * a l-V"(q) 5q-2Tp f'(q)f"(q)5q + v(t)f"(q)Sq] + 0(5q 2 ) (3.8) 

From Eq. (j3.8|) it is easy to calculate the equations of motion for quantum correction in 
the lowest order as 



J t ^ = T[r(q)} 2 ^ V " {qm2) ~ 2r P f'W"W( 5 « 2 )+vW%qm 2 )] (3-9) 

A simplified expression for the leading order quantum correction term (5q 2 ) can be es- 
timated by neglecting the higher order coupling terms in the square bracket in Eq. ()3.9|) 
and rewriting it as d{Sq 2 ) = r ^ q ^2 V"(q)(bq 2 ) dt. The overdamped deterministic classical 
motion on the other hand gives dq = — v ^T^p dt. These together yield after integration 

(5q 2 ) = A g [V'(q)] 2 (3.10) 
where A q = f^I \°p and q refers to initial position. 
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B. The stationary solution of Fokker-Planck equation in overdamped limit 



The Fokker-Planck equation Eq. ()3.5|) in the overdamped limit can be rewritten in a more 
compact form as 



dP(q,t) d 



dt 



dq T h(q) 



Do d g{ q y 



P(q,t) 



(3.11) 



T dq h(q) 

Eq. ()3.11|) is the quantum Smoluchowski equation for multiplicative noise (damping is 
corrected upto 0(1/T)). To capture the essential content of state-dependent diffusion we 
now proceed as follows. Under the stationary condition 



dP(q,t) 

dt 







(3.12) 



Eq. (j3.11J) reduces to 



D d 
T dq 



h{q) 



Pst{q) 



V\q) - Qv = 



(3.13) 



After integration of Eq. ()3.13j) we have the stationary probability distribution in the over- 
damped limit as 



P s t(q) = N- 



[9{q) 2 /h{q)\ 



exp 



i V' (q f ) 

quan V^f / 

D(q') 



dq' 



(3.14) 



with 



D(q) = frw (fi(u ) + ^ [g{qf/h{q)] 

and 



Vquan = V(q) - Q V 

and iV is the normalization constant. In the classical limit we have g{q) 2 / h(q) = 1, Qv = 
and hwo(n(cuo) + |) = k B T; the stationary probability distribution function (|3.14|) reduces 
to classical equilibrium Boltzmann distribution P^ t = N exp[—V(q)/kB r T]- Therefore the 
stationary distribution (|3.14|) is essentially a quantum generalization of Boltzmann factor for 
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state-dependent diffusion. This diffusion arises where the inhomogeneity is due to quantum 
corrections entangled with nonlinearity of the system-bath coupling. The state-dependent 
diffusion is wellknown in several classical contexts as pointed out earlier notably in Landauer 
Blow-torch effect^ 2 , where localized heating in a region along the reaction coordinate lying 
between the lower energy minimum and barrier top can raise the relative population of the 
higher minimum above what is allowed by Boltzmann factor. The example also includes 
noise-induced transport processes due to phase difference between periodic modulation of 
drift and modulation of diffusion as suggested by Biittiker 4 . It is however worth-noting that 
the inhomogeneity of D(q) as implied in Eq. 1)3.14)1 is reminiscent of some sort of quantum 
nonlocal effect. A number of points are pertinent here. First, in realizing inhomogeneity 
through D(q) we have taken care of quantum corrections to all orders. Second, the state- 
dependence is essentially due to nonlinear coupling mechanism in Hamiltonian (j2.1j) and is 
therefore model dependent and because of c-number nature of our treatment it corresponds 
to one of the classical forms of the phenomenological state-dependent diffusion term. This 
correspondence makes classical-quantum correspondence more clear. It is also worthnoting 
that the mo del- dependent nature of escape rate in classical Blow-torch effect had been 
discussed recently^. Third, for harmonic potential and linear coupling f(q) = q, we have 
Qv = and g{q) 2 /h(q) = 1 and the probability distribution function reduces to well-known 
Wigner distribution for harmonic oscillator V(q) as Pg t (q) = iVexp[— V(q)/hw (n(uJo) + |)]. 
Fourth, although model-dependent the inhomogeneous nature of the quantum diffusion must 
be thermodynamically consistent. To this end we examine in the next section the equilibrium 
condition as well as typical nonzero current situation (due to symmetry breaking under 
special condition) for a periodic potential and periodic derivative of coupling function. By 
thermodynamic consistency we imply that in absence of any external field, no directional 
component should remain after appropriate averaging over ensemble or over the period of 
space or time. 
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IV. APPLICATION: A PERIODIC POTENTIAL 



A. Solution under periodic boundary condition; Thermodynamic consistency 



In the overdamped limit the stationary current from Eq.(3.11) can be represented as, 



J 



Th(q) 



v{q) - Qv + ^d- q [m) 



PsM) 



(4.1) 



Integrating the Eq. (j4.1j) we have the expression of stationary probability distribution in 
terms of stationary current as, 



PsM) 



-Hi) 



WflKq)] 



9(0? 
h(0) 



u o Jo 



(4.2) 



where 



u o Jo 



q V'(q') - Q\ 



lo WY/Kq')} V (4 ' 3) 

We now consider a symmetric periodic potential with periodicity 2tt, i.e., V(q) = V(q + 
2tt) and periodic derivative of coupling function with the same periodicity as that of the 
potential, i.e., f'(q) = f'(q + 27r). 

Since the potential is periodic, Qy is also a periodic function because Qy = V(q) — (V(q)). 
Similarly Qf is also periodic because Qf = (/'(<?)) — f'(q), and also Q3 + 2f'(q)Qf is also 
periodic since Q 3 + 2f'(q)Q f = (lf'(q)] 2 ) - \f\q)} 2 - So from Eqs.Q and Q it is clear 
that h(q) and g(q) are also periodic functions of q with periodicity 2ir. This will be made 
more explicit when we consider a specific example in the next section. 

Now applying the periodic boundary condition on P s t(q), i-G-, Pst(q) = P s t(q + 27r) we 
have from Eq. ([4.2)1 



g(Q) 2 

h(0) 



PsM = J 7T [1 - e^] / h(q) e^dq 
^o Jo 



(4.4) 



By applying the normalization condition on stationary probability distribution which is 
given by, 



2tt 



(4.5) 
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we obtain from Eq. (j4.2j) 



p2tt p -<f>{q) [ n (C\V F2 fQ / i\ ~ 

I wfim m PJ0) - J o,l m ^ )dq '\ dq = 1 (4 - 6) 

Elimination of P s t(0) from (|4.4jl and ()4.6|) yields the expression of stationary current 
after some rearrangement 



Do 1 - 



12 JiT ^ e " 0(<?) ^ C - [1 - e«W] C e ~ m So ^Wdq 

(4.7) 

From the condition of periodicity of potential and different quantum correction terms 
it is clear that for the periodic potential and the periodic derivative of coupling function 

with same periodicity [ g (^yh(q)] anc ^ [g(g)VfeL)] are ^°^ n P er i°dic with same periodicity. This 
makes effective potential 4>(q) equal to zero; 

Therefore the numerator of the expression for current (|4.7jl reduces to zero. We thus 
conclude that there is no occurrence of current for a periodic potential and periodic deriva- 
tive of coupling with same periodicity. At the macroscopic level this confirms that there is 
no generation of perpetual motion of second kind, i.e., no violation of second law of ther- 
modynamics. Therefore the thermodynamic consistency based on symmetry considerations 
ensures the validity of the present formalism and of the overdamped quantum multiplicative 
Langevin equation. 



B. Phase induced current 

Biitikker— in 1987 had shown that in case of space-dependent friction in the overdamped 
limit a classical particle under a symmetric sinusoidal potential field and also in presence 
of a sinusoidally modulated space-dependent diffusion with same periodicity experiences a 
net drift force resulting in generation of current. This current is basically due to the phase 
difference between the symmetric periodic potential and the space dependent diffusion. The 
current does vanish when the phase difference is either zero or integral multiple of tt. van 
Kampen= in a later work also ended up with similar kind of conclusion for a system with 
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space dependent temperature under the overdamped condition. The result of van Kampen 
is a re-examination of the earlier observation due to Landaueri^. 

We now explore, in the spirit of Buttiker but in the present quantum-mechanical context 
a system with a symmetric periodic potential, i.e., V(q) = V(— q), and a periodic derivative 
of coupling function with same periodicity but with a phase difference between them leading 
to a net directed motion or current. This is because of the fact that the phase bias gives a 
tilt to the effective potential <fr(q), (if <f)(q) is plotted as a function of q), which makes the 
transition between left to right and right to left unequal. The phase difference therefore 
breaks the detailed balance of the system. We show that when the phase difference is zero 
or integral multiple of 7r the quantum current vanishes. 

We proceed with a sinusoidal symmetric potential of the form 

V(q) = V [1 + cos(g + 9)} (4.9) 

where Vq can be taken as the barrier height and 9 is phase factor which can be controlled 
externally. 

The derivative of the chosen coupling function, f(q) = (q + 0! sing), is 

f\q) = 1 + acosg (4.10) 
where a is a modulation parameter. 

Given the form of potential V(q) and coupling function f(q) we obtain from Eq. ([3.100 
the second order quantum correction in the overdamped limit as 

(5q 2 ) = - A q V 2 sin 2 (q + 9) (4.11) 
Therefore the correction to the potential in the leading order from Eq. ([2.2fi|) is given by 

Q v = -±A q V 3 sm 3 (q + 9) (4.12) 

The quantum corrections Qf and Q 3 in the same order can be estimated using Eqs. ([2.29[) 
and flODj) : 



Qf = -- A ? q cosg sin 2 (g + 
Q 3 = A q a 2 V 2 sin 2 g sin 2 (q + 9) 



(4.13) 
(4.14) 
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Furthermore from Eqs. (j3.3j) and ()3.4|) we calculate h(q) and g(q), respectively, using 
Eqs.flTinj), KT^i and KWi to obtain 



h(q) = (1 + a cosg) 2 — A q a Vq cosg sin 2 (g + 9) (1 + a cosg) 

+ A q a 2 V 2 sin 2 q sin 2 (g + 9) (4.15) 

g(q) = 1 + a cosg — - A q a V Q 2 cosg sin 2 (g + 9) (4.16) 

With these preliminaries we now calculate the current given by Eq. (j4.7|) explicitly using 
IjOj) . (f4~T2l. (|4~THll and (OHl) . A key quantity for this analysis is 0(g) [ (|Ojl ] which involves 
the phase In units of h = &b = 1 we set the parameter values (5q 2 )o = 1/2, minimum 
uncertainty value, A q = 0.5, Vq = 1.0, ujq = 1.0, a = 1.0, T = 1.0 and T = 1.0. The variation 
of current J as a function of phase 9 is exhibited in Fig.l. It is interesting to observe that for 
6 7^ 0, mr where n = ±1, ±2..., phase induces a current which is a periodic function of the 
phase difference between modulations of potential and diffusion. The origin of the current 
can be traced to the inhomogeneous diffusion of a quantum particle in contrast to classical 
one proposed by Biittiker. That this current depends on the amplitude of modulation a of 
the diffusion is shown in Fig.2 for 9 = 3.6, Vq = 1.0, ujq = 1.0 for several values of ksT. 
For a = the current vanishes and we have only linear coupling and the phase bias has no 
relevance in such situation. We observe for a fixed temperature a sharp increase in current 
beyond a moderate value of a and that for an optimum temperature the current is maximum 
for a given strength of modulation. Fig. 3 illustrates the variation of current (J) as a function 
of temperature for the phase 9 = 3.6 and for several values of strength of modulation of 
diffusion a for Vq = 1.0, uoq = 1.0. One observes that even at T = 0, the vacuum field of the 
heat bath induces a finite current. As the temperature increases the current decreases after 
an initial increase and reaching a maximum. 

We thus observe that nonlinear system-bath coupling may give rise to state-dependent 
noise and diffusion in a quantum system. For a periodic potential and for a periodic deriva- 
tive of coupling function, with same periodicity, the state-dependent noise may lead to 
symmetry breaking in presence of a phase bias. This generates a directional flow which 
vanishes in absence of the bias and can be optimized by application of suitable strength of 
modulation and temperature. 
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V. CONCLUSION 



In this paper we have developed a theory of diffusion of a quantum particle in inho- 
mogeneous media. The approach is based on the system-reservoir model with a nonlinear 
coupling. We derive the quantum Langevin equation with a multiplicative noise and a non- 
linear dissipation in the Markovian limit, which is coupled to a set of quantum correction 
equations developed order by order. A systematic expansion of the relevant variable in pow- 
ers of inverse of the dissipation constant and use of large friction limit lead to a quantum 
Smoluchowski equation for state-dependent diffusion. It is apparent that the state depen- 
dence owes its origin to nonlinear coupling between the system and bath degrees of freedom 
and the corresponding generalization of Boltzmann factor for the steady state has been 
shown to be thermodynamically consistent. We have applied the formalism to the problem 
of diffusion of a quantum particle in a periodic potential where the derivative of coupling 
function is periodic with same periodicity. We have shown that a phase difference between 
these two spatially periodic modulations may give rise to a directed quantum current. This 
current vanishes in the classical limit and is a consequence of state-dependent diffusion where 
nonlocality in the effective potential is essentially a quantum effect. 

The Brownian motion of a quantum particle in inhomogeneous media is an active area 
of contemporary research in stochastic energetics. We hope that our formalism of quantum 
Langevin equation with multiplicative noise and the associated Smoluchowski equation for 
state-dependent diffusion as a description of the stochastic processes may be useful for 
other similar issues, particularly for calculations of quantum decay rate of metastable state, 
energy dissipation and so on. The extension of the theory to non-Markovian and weak 
friction regime is also worth-pursuing. We hope to address these issues elsewhere. 
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Figure Captions 



Fig.l: Variation of current J as a function of phase 6 between to 2n for T = 1.0, 
lu = 1.0, a = 1.0, V = 1.0 and T = 1.0 (units arbitrary). 

Fig. 2: Plot of current J vs. strength of modulation a for 9 = 3.6, V = 1.0, u = 1.0 
and for T = 0.1 (dashed-dotted line); T = 0.5 (dotted line); T = 1.0 (dashed line); T = 3.0 
(solid line) (units arbitrary). 

Fig. 3: Variation of current J as a function of temperature T for 6 = 3.6, Vq = 1.0, 
uo = 1.0 and for strength of modulation a = 0.5 (dashed-dotted line); a = 0.7 (dotted line); 
a = 0.9 (dashed line); a = 1.0 (solid line) (units arbitrary). 
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